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Abstract 

We show that the vacuum structure of a generic multi-Higgs-doublet model shares several impor- 
tant features with the vacuum structure of the two and three Higgs-doublet model. In particular, 
one can still define the usual charge breaking, spontaneous CP breaking and normal (charge and 
CP preserving) stationary points. We analyse the possibility of charge or spontaneous CP breaking, 
by studying the relative depth of the potential in each of the possible stationary points. 

PACS numbers: ll.30.Er, 12.60.Fr, 14.80. Cp, 11.30.Ly 
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I. INTRODUCTION 



Most features of the Standard Model (SM) of electroweak interactions have been probed 
to a very high precision. Still, the Higgs sector remains largely untested and new physics is 
certainly possible. In particular, one might have more than one Higgs, as required for ex- 
ample by supersymmetry. Multi- Higgs models are also appealing for a variety of theoretical 
reasons related to CP violation: i) if the Higgs potential conserves CP, this symmetry could 
be spontaneously broken by the vacuum; ii) if there are three or more Higgs doublets, then 
there might be CP violation in the mixing matrix of the charged Higgs; and iii) the new 
sources of CP violation have the potential to explain bariogenesis. 

A drawback of these models is that they involve a large number of parameters. For 
example, the scalar potential of the most general two-Higgs-doublet model (2HDM) involves 
14 real parameters, while the potential which explicitly preserves CP involves 10 parameters 
(these may be reduced to 11 and 9, respectively, through suitable basis choices). Given 
the large number of parameters present in these models, a variety of methods have been 
developed in order to restrict the parameter space, some related to the vacuum structure of 
the scalar potential. Recently, some interesting features of the vacuum structure have been 
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obtained for the particular case of the 2HDM [1]. Namely, it was shown that, whenever a 
normal - charge and CP conserving - minimum exists in the 2HDM, the global minimum of 
that potential is the normal one. Moreover, it was shown that the depth of the potential 
at a stationary point that breaks charge or CP, relative to the normal minimum, is related 
with the squared mass matrix of the charged or pseudoscalar Higgs (evaluated at the normal 
minimum), respectively. Recent work on these subjects may be found in |2]. In this work we 
will analyse how these conclusions may, or may not, be generalised to the case of a potential 
with N Higgs doublets. 

The paper is organized as follows. In section ^ we introduce our notation and prove one 
of our main results: that the vacuum structure of a generic multi-Higgs-doublet model may 
be reduced to vacua involving only two or three doublets. This is accomplished through a 
series of basis transformations, for which the potential is invariant, even if its parameters are 
not. We discover that the study of possible charge breaking (CB) vacua is more easily done 
in a basis where only three out of the N doublets have non-vanishing vacuum expectation 
values (vevs). For CP violation, the appropriate basis is even simpler: only two doublets are 
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non-zero. In section ITTTl we compute the values of the potential at the charge breaking and 
normal vacua and compare their values. We study whether it is possible to obtain charge 
breaking minima deeper than a normal minimum. In section IIVI we repeat this procedure, 
but now for CP breaking vacua. We present our conclusions in section |3 Appendix 1X1 
provides a basis and gauge independent definition of the CB vacuum, while appendix ^ 
contains a specific example of a three Higgs-doublet model for which the CB vacuum lies 
below the normal vacuum. 



II. THE SCALAR SECTOR OF A GENERIC N-HIGGS-DOUBLET MODEL 



A. The scalar potential 



ow closely the notation of Refs. yl yi, where more details may be 

nnnnn 

p, la, m la l2| • Let us consider a SU(2) 011(1) gauge theory with N Higgs- 
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found - see also 
doublets with the same hypercharge y = 1/2, denoted by 





hi 

















where Vi are their vacuum expectation values (vevs), and i runs from 1 to N. In all that 
follows, we will use the standard definition for the electric charge: Q = + Y , meaning 
that all vevs in the lower components of the doublets are electrically neutral. With this 
definition, a vacuum with all upper components of the vevs equal to zero, z/" = 0, does not 
break the charge symmetry. 

The scalar potential may be written as 

Vh = + A.,,feK$l$,)($I$0, (2) 

where Hermit icity implies 

* 

Kj,ki = ^ki,ij = ^*ji,ik- (3) 

Under a unitary basis transformation of the Higgs fields, their kinetic terms remain the same 
but the coefficients fiij and Xij^ki are transformed in such a way that the potential remains 
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invariant. Using Eq. (^J, the scalar potential becomes 

t 



+z/.' Uiij + 2\ij^ki vlvi + ifl fiij + 2\ij^ki vlvi Vj 
{^yj) + 2\ijM (^l^j) {j^Ij^i) + 2\a,kj {J^lvj) 
+\ijM (^l^^i) {vi^i) + \ijM {T^lvj) i^lvi) 
+2Xij,ki ivlvj) ivli^i) + 2Xij,ki ivlvj) (J^lvi) 
+^ijM ivlVj) ivWi)- 



(4) 



Requiring in Eq. (0} that the linear terms in (fi vanish, gives us the stationarity conditions 



(forz = I,-- - ,N). 



(5) 



ki ^k^i 

Multiplying by uj leads to 

f^iji^l^j) = -2Aij,H (z/|z/j) (z/Jz/;). (6) 

The value of the potential at a stationary point is found from Eq. Q by setting all (pi = 0. 
Using Eq. ©, this may be written in the following three forms 



-KiM ivlvi) 



(8) 
(9) 



B. A simple basis to study charge breaking in the NHDM 



After spontaneous symmetry breaking the Higgs fields acquire the vevs 



(10) 



An analysis of the potential with such a complicated vev structure would be too difficult 
to perform. We will now show how, using the freedom to choose a basis for the Higgs 
doublets, one manages to simplify immensely this study. We start by performing a unitary 
transformation on the last two Higgs fields according to 



With this transformation, the vevs of the last two fields become 

\,,u |2 I I , ,M 12 



^^^^ /\7/u 12 _|_ 12 I u d I u d ^ ' ^''"^■^ 

respectively. We have thus succeeded in removing the upper component of the vev of the 
last Higgs field. Moreover, the upper component of became real and positive. We 

can continue with similar transformations, applied to successive pairs of Higgs fields, until 
the corresponding vevs become 



\ / \ / 



(13) 



Notice that z/" and vf in Eq. are not the same as in Eq. but rather the values 

obtained after the successive transformations of the type shown in Eq. (fTT|) . Similarly, we 
keep the notation for the fields, although they have been transformed through a series of 
basis changes. At the end of the process outlined, is real and positive. 

We may now repeat the exercise with the lower components. Indeed, through the trans- 
formation 



we can change the vevs of the last two Higgs fields in Eq. (fT!?|l into 





T,d 12 I I , ,d 12 



^'n) = \ , (15) 



respectively. This eliminates the lower component of (^'^v) ^^"^ makes the lower component 
of real and positive. We may continue with the other down components, until we 

reach the following vev structure 







(16) 



Notice that li^al, which is real and positive due to Eq. (fT^ . cannot be removed without 
implying the appearance of an upper component on the second vev. We have thus reached 
a simple but remarkable result. Indeed, althou gh t here are many parameters involved in 
the general N-Higgs doublet model, its vacuum 10| structure can be brought into a much 
simpler form, through a suitable basis choice. 

If, after all these basis transformations, we are left with a vev structure for which ^ 0, 
then the vacuum breaks electric charge. As in the 2HDM, we may now utilize the gauge 
freedom in order to bring the vevs into the final form |1]J 



Vcl / \ Vc2e"''' 





(17) 



where 5c is a phase, while a, fd, and Vc2 are positive real numbers. This, then, is the 
simplest form one can find for a CB vacuum. 

However, we are interested in comparing the value of the potential at the CB vacuum 
with its value at the normal vacuum. Therefore, we must find out what is the form of the 
most general normal vacuum, in the basis in which Eq. (fT7j) is written. Clearly, given our 
definition of electric charge, it will have all = 0. A generic charge-preserving vacuum, 
then, will have the form 





\ / 



'N-l / \ ''Af 



(18) 



We emphasise that these vf are not the same as those appearing in Eq. (fTUj) . However, 
we can now apply the same method we used previously to bring the normal vacuum to a 
more manageable form. Through a transformation analogous to that of Eq. ()14|) . we can 
set the last doublet to zero. Notice that this basis change does not involve the first two 
doublets, so the charge-breaking vacuum structure, Eq. (fTTj) . remains unaffected. Successive 
basis transformations may be applied that do not change Eq. (fTTj) but set to zero the lower 
component vevs of Eq. (fTHj) . until one is left with the final normal vacuum structure. 



\ / 






(19) 



If we try to perform another basis change to set to zero, we will destroy the simple 
form for the charge-breaking vevs, Eq (|T7jl. We denote by the "B-basis" the basis where the 



charge breaking vevs have the simple form of Eq. ()17|) and the normal vacuum vevs, the form 
given by Eq. ()19|) . The B-basis is appropriate to study the possibility of charge breaking 
vacua. This result shows that the study of charge breaking for an N-doublet potential is 
reduced to the analysis of the three-doublet situation. 



C. A simple basis to study CP breaking in the NHDM with explicit CP conser- 
vation 

We are also interested in the possibility of CP being spontaneously broken in N-Higgs 
doublets models. In this case we cannot simply choose the most general NHDM potential 
- we must make sure that that potential does not break CP explicitly. In appendix A of 
0], Gunion and Haber invoke CPT and = 1 (where T is the time- reversal operator) to 
show that: "The Higgs potential is explicitly CP-conserving if and only if a basis exists in 
which all Higgs potential parameters are real" . We therefore consider one such basis for our 
NHDM potential: all of its parameters are real and it explicitly preserves CP. 

Given our definition of electric charge, the most general charge-preserving vacuum (CP 
violating or not) will be of the same form as Eq. ()18|). Our starting point, however, is not the 
basis in which we wrote Eq. ()18|) . but rather a generic basis for which the parameters of the 
potential are all real. It is now convenient to ensure that all basis changes do not introduce 
complex parameters in the potential. This means that we are restricted to orthogonal basis 
transformations. Even with this restriction we are still able to simplify immensely the study 
of the NHDM potential. This is accomplished through two series of steps: 

1. We start with an orthogonal transformation on the last two Higgs fields according to 

•^-^-i \ ^ 1 ( -Im -Im (i^i) \ I 

<^'n I + W «) \ -Im «) Im J ) \ <^n 

(20) 

This eliminates the imaginary part of the vev of the last Higgs field. We can con- 
tinue with similar transformations, applied to successive pairs of Higgs fields, until the 
corresponding vevs reach a structure of the type 










(21) 
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Notice that the vf in Eq. ()2ip are not the initial ones, but rather the values obtained 
after the successive transformations of the type shown in Eq. ((201) • Similarly, we keep 
the notation for the fields, although they have been transformed through a series of 
basis changes. After these steps Im [vf) < 0. 

2. We continue with the orthogonal transformation on the last two Higgs fields 

•^-'^-i \ 1 Re Re 




(22) 

^'n ) ^Re^ «_0 + Re^ «) V "^^^ K) ' ' ' 

This eliminates the vev ($'^), simultaneously making the lower component of 

real and positive. We can continue with similar transformations, applied to successive 

pairs of Higgs fields, until the corresponding vevs reach a structure of the type 

(23) 








where Re (z/g) > 0. 

Hence, when the scalar potential conserves CP, it is possible to choose a basis in which 
only the first two doublets have vevs, while keeping all parameters in the potential real. At 
this point, we distinguish two physically distinct scenarios. If lm{uf) =0, then the vacuum 
is a normal one and preserves CP; if not, it spontaneously breaks that symmetry - we call 
it a CP violating (CPV) vacuum. 

Let us then suppose we started with a normal vacuum, and that we employed the basis 
transformations described above until the vacuum structure was reduced to Eq. ()2H|1 . with 
real vevs uf = vi and z/g = f2- Because the remaining vevs are real, we can perform a final 
basis transformation on the first two fields 



bringing their vevs into the form 



$2 / Vvf + vl 



1 Vl V2] 


II 




\ -V2 Vl j 




^4. J 



(24) 





^0^ 






ni = 1 






: 




iv) 







(25) 



where v = \/vl + v^, and all remaining doublets are zero. This is known as the "Higgs 
basis" for the normal vacuum in the 2HDM 
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However, we are interested in comparing the value of the potential at the normal vacuum 
with its value at a CPV vacuum. Therefore, we must find out what is the form of the most 
general CPV vacuum, in the basis in which Eq. (|25p is written and with the definition of 
electric charge we have adopted. It will be of the form of Eq. (jl8|) . with new vevs uf, 



(26) 

























\Z2 



Because the normal vacuum has been reduced to the form of Eq. p5|l . where only the first 
doublet is different from zero, we can apply the steps 1 and 2 detailed above for this new 
vacuum and again reduce the CPV vacuum to the form 

(27) 

with all remaining doublets zero. In this equation, Zi and Z2 are complex numbers. The 
sequence of steps that led us from Eq. ()26|) to Eq. (jTTj) did not change the normal vacuum of 
Eq. (j25|) . because none of those steps involved the first doublet. We cannot further remove 
the imaginary or real parts of Z2 because that operation would involve the first doublet and, 
thus, take us away from a basis in which the form of the normal vacuum of Eq. remains 
valid. The basis for which the normal and the CPV stationary points have the simple form 
given by Eqs. ^ and ^ will be called the "S-basis" . 

The S-basis is very useful because, when using it for the normal minimum, the Goldstone 
bosons are isolated as the components of oji, while the other iOi (i = 2, ■ ■ ■ , N) contain other 
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charged and neutral scalars fields p|,y, yj. Indeed, in the S-basis 

ipi = { I , (for z = 2,--- ,iV), (29) 

iR, + tI,)/V2 J 

where and are the Goldstone bosons (which, in the unitary gauge, become the 
longitudinal components of the and of the Z^); couples to fermions proportionally 
to their masses (in the fermion mass basis); and H^, Ri, and /, {i = 2, ■ ■ ■ , A^) are the 
charged and neutral scalars fields. Notice that these are not the physical particles; those 
will be obtained by diagonalizing the squared mass matrix of the charged Higgs, and the 



squared mass matrix of the neutral Higgs (including H'^, Ri and li). These important 
properties will become obvious below. 



D. The mass terms at the normal vacua 

We now wish to study the quadratic terms in Eq. when the vevs are taken to coincide 
with those at a normal vacuum. Since the basis transformations do not mix the upper and 
lower components, the normal vacua have z/" = for all i, in any basis (as long as no gauge 
transformations are made). As a result, the quadratic terms of Eq. (0)) evaluated at a normal 
stationary point may be written, in any basis, as 

iM% = {Ml)l + {Ml)l Re {^f) Re (^?) + {M^)l Im [^f) Im (^^^) 

+ {Mh)l Re (^f) Im (y.?) + (M|^)^. Im (^^) Re (^?) , (30) 

where we identify 

{Ml)l = /i,,+2A,,,H//fVf, (31) 

{Ml)l = Re [/X,, + uf + 2A,,,,,z/,^ uf* + 2\u,A A , (32) 

{M^)l = Re [fi,, + 2A,,- uf + 2A,,,,,z/,^ uf* - 2A,,,,,z/,^ z/f] , (33) 

{Mh)l = -Im K- + 2A.,,H^f + 2A,,,,,z/f z/f * - 2A,fc,,,z/f z.f] , (34) 

(M,2^)^. = Im [ft,, + 2X,,,,iut* uf + 2\,^,,ut uf* + 2A,fc,,,z/,^ uf] , (35) 

and the superscript N indicates that these mass matrices have been evaluated at the normal 
vacuum. 

Using Eqs. 0, one can show that the matrix (M^)'^ is hermitian, while the real matrices 
(M^)^ and [M'j)^ are symmetric. The remaining two matrices are real and related by 
(M|^)^. = {Mj^)^.. This implies that the 2N x 2N matrix. 



(M|)N (M2,)N 

2 \N ( T\/T2\n 



(36) 

{M^r {Mr ' 

is symmetric. Moreover, the matrix [M"^^ behaves like a second rank tensor under a basis 
transformation of the Higgs fields, but the other matrices do not. 

As we mentioned, these expressions are valid for normal vacua in any basis. For the 
B-basis, where only {vf , 1/3 ) ^f} ^i-re different from zero, the indices {k , /} in the Eqs. ()31|) - 
run only from 1 to 3. In the S-basis, where only the first doublet has a non-zero vev, 
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the mass matrices are simplified considerably. In what follows the vevs, the parameters /ijj, 
and Xij^ki are all written in the S-basis; it is important to understand that changing the basis 
would change the vevs, but also the parameters fj,ij and Xij^ki [S]- We find, 

{Ml)l = fi,, + 2v'X,,,u, (37) 

(M2)[^. = Re [/i,, + 2v^ (Ai,- n + Aii.ii + X^l,Jl)] , (38) 

(M2)[^. = Re [/i,, + 2v^ (Ai,- n + Ki^ij - Kiji)] , (39) 

{mIj)^^ = -Im [fiij + 2v^ (A,,- 11 + Aii,ij - Aiiji)] , (40) 

ij 



{M]ji)l = Im [fi,j + 2v^ (Aij, 11 + Xii^ij + Aiiji)] . (41) 



Using the parametrization of the normal stationary point in the S-basis, shown in Eq. ()25|). 
on the stationarity conditions of Eq. (0), we find 

/i,i + 2v^X^i,n = 0. (42) 

But this coincides with the definition of (M|)^^, in Eq. ()37p and, since this is a hermitian 
matrix, we conclude that the first row and the first column of {M]_)^ have zero in every 
entry, 

{Ml)l = = {Ml)l (for^ = l,---,iV). (43) 

This shows that, indeed, in the S-basis coincides with the charged Goldstone boson, in 
accordance with Eq. (j2Hl)- 

Also, using Eq. (fH^ with j = I and Eq. 



{M^)l = Re {fia + 2v^Xa,u) + 2v''Re (A,i,ii - A,i,ii) = 0. (44) 

Since (Mj)^ is symmetric, we find 

{M])l = = {M])l (for^ = l,---,iV). (45) 

To simplify, let us now consider for a moment the case in which all fiij and all Xij^ki are 
real. As we explained above, a CP-conserving NHDM potential falls under this category. In 
that case, (M^^)'^ = = (M|^)^, and the matrix in Eq. becomes block diagonal. In 
addition, (M^)^ and are the squared mass matrices of the scalars and pseudoscalars, 

respectively. Thus, Eq. ()45|) shows that in the S-basis lm{(ff) coincides with the neutral, 
pseudoscalar Goldstone boson, in accordance with Eq. 
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III. THE CHARGE BREAKING VERSUS THE NORMAL STATIONARY 
POINTS 



Throughout this section we will work in the B-basis, although it will be obvious that our 
final results hold in any basis. We assume that both the N and CB stationary points exist. 
We will now compute the values of the potential at each of those stationary points. To that 
effect, we first recall the results of section |nl where we showed that, in the B-basis, the vev 
structure of both stationary points is given by 



Charge breaking vacuum (CB), 





f a \ 




( ^ 
















\ I'd 1 















ni = 






I -J 












ns = , Normal vacuum (N), (46) 



with all remaining doublets having vevs equal to zero. The parameter Sc is a phase, while 
a, Vci and vs are real, positive numbers. In general, ui and z/2 are complex. From Eq. (jSJ, 
we obtain the stationarity conditions for the CB stationary point in the B-basis, 

(^/iii + 2Xa,ki 4q) a = (47) 

(^fiii + 2Aii,H clci^ Vci + (^fii2 + 2Ai2,fc« 4c/ j Vc2e'^'' = , (48) 

for i = 1 ■ ■ ■ , N and k,l = 1,2. Since we assume that the CB stationary point exists, 
a ^ and its coefficient in Eq. (|47p must equal zero. From Eq. (|48|) . then, the coefficient of 
Vc2e^^'' is also zero. As a result, the stationarity conditions at the CB stationary point may 
be written as 

fiij, + 2\j,^ki 4q = (i = 1 ■ ■ ■ , AT ; ji = 1, 2). (49) 

Let us now contract the indices {i , ji} with n\nj^. This gives 

(^/ijj, + 2Ayi,fc, 4q j n\nj^ =0 (z = 1, 2, 3 ; ji = 1, 2), (50) 

and from here it is trivial to obtain 

/iy n\nj + 2Aij- M n\njclci - (^fi^ + 2Ai3,fc/ 4^;) nln^ = (z, j = 1 ■ ■ ■ , 3). (51) 
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Notice the appearance of the term ^ijn\nj which, according to Eq. (jH}, equals twice the 
value of of the potential at the normal stationary point, . 

Now, from Eq. (j3ip . the mass matrix for the charged scalars at the N vacuum in the 
B-basis is given by 

{Ml)l^ = iJki + 2Xki,^jnlnj . (52) 
Contracting the indices {k , 1} with c|.q we obtain 

(^±)m = /ifc/ clci + 2Xki,ij nluj eld 

= + 2\i,^kin\n,c\ci , (53) 

where we have used Eq. dH} to identify /ifc^ c|.q as twice the value of the potential at the CB 
stationary point and the symmetries of the A coefficients from Eq. Comparing Eqs. (|51|) 
and we can subtract them to find 

VS'' -VS = \ {Ml)l c\c, -\v^ (/i.3 + 2A,3,fc^4Q) . (54) 

This is our main result regarding the possibility of charge breaking in the NHDM. Although 
obtained in the B-basis, it is very simple to rewrite Eq. ()54|) in a basis invariant form. At 
this point it is important to recall the results obtained in reference |Q for CB in the case of 
the 2HDM. In the notation of this paper, the conclusions therein reached are written as 

y^"" = \ {Ml% c\c, . (55) 

When the normal stationary point is a minimum, the matrix (M|)^ has, besides the Gold- 
stone bosons, only positive eigenvalues, and it is very easy to prove |l| that one obtains 
— > 0. Hence, if a normal minimum exists, the CB stationary point is always 
above it. No possibility of tunneling from the normal minimum to a deeper CB stationary 
point exists in the 2HDM. 

The similarity with the NHDM case is clear, but the difference of the potential depths 
now contains an extra term, proportional to ^3. Let us consider that the normal vacuum in 
he NHDM is indeed a minimum. Then, as before, the term (M^.)^. c|cj is strictly positive 
1^ . However, there is no a priori reason for the second term in the right-hand side of 
Eq. to be positive. And in fact, depending on the values of the parameters /i and 
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A, it may well be negative, so much so that it overwhelms the positive contributions from 

As an example of this possibility, we undertook a study of CB in the 3HDM for generic 
values of the parameters of the potential. For simplicity we considered the 3HDM potential 
without explicit CP violation. Our conclusions follows: 

1. As in the case of the 2HDM, it is certainly possible to find combinations of {/i , A} for 
which there are normal minima with a CB stationary point located above them. 

2. However, unlike the 2HDM situation, we have found combinations of {/i , A} for which 
both the normal and Charge Breaking stationary points are minima, hut verify < 

yN 

In Appendix IHl we give a set of numerical values of {fj, , A} corresponding to this situation. 
In fact we obtain, from such parameter values, 

Vg^ = -2.6678 X 10^ GeV^ < = -2.2792 x 10^ GeV^ . (56) 

A numerical minimization of the potential found no value below V§^. 

To ensure that both CB and N are minima, we calculated the scalar squared mass matrices 
at both stationary points. Other than the expected zero eigenvalues (3 for the N minimum, 
4 for the CB one) all the others are positive. 

In conclusion, the study of charge breaking vacua in the NHDM reduces itself to the 
study of a 3IIDM potential. And, for this one - and unlike the 2HDM case - there is the 
possibility of CB minima which are deeper than a normal minimum. 



IV. THE CPV VERSUS THE NORMAL STATIONARY POINTS 

Let us consider a Higgs potential with explicit CP conservation and no CB stationary 
points. As we showed in section Inl it is possible, through a series of orthogonal transforma- 
tions that preserve fi and A as real (though changing its values), to reach what we called the 
S-basis, where a normal (CB and CP preserving) and CPV (CB conserving, CP violating) 



14 



stationary points have vevs given by 

si = I I , S2 = I I , CP- violating vacuum (CPV), 









:;) 










, ^2 = 


(:) 







ni = , ^2 = , Normal vacuum (N) (57) 

VI \V 

where v is real and at least one of {zi , Z2} is complex. Unlike the CB case, we are able 
to reduce the study of CP violation in the NHDM to the analysis of only two doublets. 
Throughout this section we will work in the S-basis. We now assume that both the N and 
CPV stationary points exist. Eq. (j42p shows the stationarity conditions of Eq. (0) applied to 
the normal stationary point and written in the S-basis. Similarly, using the parametrization 
of the CPV vevs in the S-basis, shown in Eq. (fTTf) . on the stationarity conditions of Eq. 
we find 

(^fiii + 2\ii^ki s\si^ Zl + {^^ii2 + 2Ai2,fc« s\si^ Z2 = 0. (58) 
Specifying for i = 1 and rearranging the terms, we obtain 

Zl /ill + Z2 /il2 = -2\ll,kl slsi Zl - 2Ai2,H slsi Z2, (59) 

This can be viewed as a system of one complex (two real) equations in the two real unknowns 
flu, and /ii2. The solutions are easily obtained. One finds that 

~ 2^^^"^ ~ ^ll.llklP + ^11,21 {zlz2 + Z2Z1) + (Aii^22 ~ ^12,12 + -^12,21) 1^2^ 

= Aii,ii|2;i|^ + Ai2,ii2i^2 + A21, 11^2^1 + (A22,ii + A2i,i2 — A2i,2i) 1^2!^ 

= {^ij,U + ~ -^iljl) ^l^j- (60) 



In addition, the stationarity condition at the normal minimum, Eq. (j^zj). yields 

- ^yUii = t'^Aii^ii = Xij^iinlrij 

= (Ay, 11 + Aii.ij - Aiiji) nluj. (61) 

We conclude that 

(Aij,ii + Aji^ij — Ajiji) s\sj = (Aij,ii + Aii,ij — Ajiji) nluj. (62) 

We are now ready to calculate the difference between the value of the scalar potential at 
the CPV stationary point and the value of the scalar potential at the N stationary point. We 
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start from the definition of the pseudoscalar mass matrix (Mf)^ in Eq. ()39p and multiply 
it, respectively, by slsj and nlrij, to find 

Subtracting both lines we find 

VS'^-V- = l{M^)lsls,. (64) 

In obtaining this result we have used Eq. (jHSl), and we noticed that, according to Eq. (jHl), 

T/CPV _ 1,,.. st„. 

^i? = ^/^ii^I^i- (65) 
Furthermore, we used the fact that Eqs. and imply that 

{M])lnln, = 0. (66) 

Eq. (jMj) is the generalization of the results obtained in ref. |1| for CP violation in the 2HDM. 

It can be shown, using the general definition of (Mf)^ in Eq. (jH^ . that Eq. (jMj) is 
invariant under orthogonal basis transformations, so that in Eq. we can actually consider 
the indices {i , j} going from 1 to N. For simplicity, we evaluate it in the S-basis. As before, 
when N is a minimum, we will have 

Vr-Vj^ = l{M^)l\z,\^>0. (67) 

The result is strictly positive because the only zero eigenvalue of the matrix (Mf)^ is in 
the first line/row; the remaining sub-matrix is definite positive. This implies that all of the 
elements of its diagonal - such as {Mf)^2 " are positive. This is another advantage of utilizing 
the Higgs basis. 

Eq. (jnZj) generalizes the results obtained in Ref. [Jj for the particular case of = 2: 
whenever a normal minimum exists it is certainly deeper than any CPV stationary point. 

Notice that one can obtain a CPV stationary point which is deeper than a normal sta- 
tionary point N. This occurs for parameters such that {Mf)^^ < 0- However, in that case N 
is not a minimum (although it is a stationary point). 
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V. CONCLUSIONS 



We have studied the vacuum structure of the most general N-Higgs-doublet model. We 
have shown that, in order to compare the depth of the potential at a normal minimum with 
its depth at a CB stationary point, a basis may be chosen such that the vacuum structure 
mimics that of the 3HDM. Similarly, in order to compare the depth of the potential at a 
normal minimum with its depth at a CPV stationary point, a basis may be chosen such that 
the vacuum structure mimics that of the 2HDM. 

This great simplification allowed us to generalize the results of Q], showing that, whenever 
a normal minimum exists, it is certainly deeper than any CPV stationary point. However, 
we found one remarkable difference regarding CB: whereas in the 2HDM it is impossible to 
find CB minima below normal ones, that does not happen for the NHDM, with > 3. This 
raises the possibility of finding charge breaking bounds j3| for these potentials, which might 
improve their predictive power. Notice, however, that if the parameters of the potential are 
such that at the N minimum (in the B-basis) one has = 0, one recovers the 2HDM result 
for the NHDM potential: if such a normal minimum exists, it is certainly deeper than the 
CB one. This can be used as a sufficient condition to prevent CB from occurring in the 
NHDM. 

It is interesting to note that Eq. shows that the difference between the value of the 
potential at the CB stationary point and the value of the potential at the normal stationary 
point is related to the charged Higgs squared mass matrix. That relation is perfect for 
the 2HDM, but "spoiled" by the v^, terms in Eq. (|54|) for the NHDM. Similarly, when the 
potential conserves CP, Eq. (j64j) shows that the difference between the value of the potential 
at the CPV stationary point and the value of the potential at the normal stationary point 
is related to the pseudoscalar squared mass matrix. Thus, the depth of a potential at a 
stationary point that breaks a given symmetry, relative to the normal minimum, is related 
to the squared mass matrix of the scalar particles directly linked with that symmetry. 
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APPENDIX A: CHARGE BREAKING - THE KINETIC TERMS 



In this appendix wc present a basis and gauge independent definition of the charge break- 
ing vacua. This could be done by looking at the Goldstone bosons in the scalar mass matrix, 
but it is easier to look at the mass of the photon instead. The kinetic terms for the scalar 
fields are 



(Al) 



where W^^ (a = 1, 2, 3) and 5^ are the SU(2)i and U(l)y gauge bosons, respectively, and 
(a = 1, 2, 3) are the Pauli matrices. After spontaneous symmetry breaking we obtain mass 
terms for the gauge fields given by 



g'B, + gW^ V2gW;: 



vf 



where 



V2 ■ 



After some reorganization, the result is proportional to 



I W"^ WP ^B, ] M?^ 



( \ 



where 
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(A2) 



(A3) 



(A4) 



(A5) 



19 



for each doublet $ j = 1, ■ ■ ■ A^). Summing over all doublets and introducing the complex 
vectors 



the mass matrix may be \ 

N 
i=l 

with the notation 

N 

{zu.Zdr = Zd.z^ = Y,{^fTK- (A8) 

i=l 

The determinant becomes 

det = 4 + |^^|2)2 [|^^|2|^^|2 _ |^^_^j2] _ (A9) 

This allows us to define the charge-breaking (CB) vacuum in a completely basis and gauge 
independent fashion. Indeed, in order for the vacuum to conserve charge, i.e., to conserve 
U(l)cimg; we need to have a massless photon. But that implies that the determinant in 
Eq. ()A9|) must vanish. Any combination of vevs {z/f , vf} for which this does not occur 
is therefore a CB stationary point. As is easily seen from Eq. ()A9|) . we cannot get a CB 
vacuum with only one Higgs doublet, a well known result. 

APPENDIX B: A THREE HIGGS DOUBLET POTENTIAL WITHOUT EX- 
PLICIT CP VIOLATION 

In this appendix we give a specific example of a 3HDM potential without explicit CP 
breaking for which one finds a CB minimum deeper than a normal one. Since no explicit CP 
breaking occurs, we work in a basis where all {/x , A} are real. The values of the parameters 
are given below. 

All remaining parameters are obtained from these using the symmetries of , A} ex- 
pressed in Eq. (jSJ. As mentioned this set of parameters gives us a normal minimum and a 
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Zd= [vf^ 



,d ,,d 



V, 
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V 



} 



N \ 1 



(A6) 



:itten as 



\Zu^ + 




^Zu.Zd 



\Zu\ + \Zd\ 









\/2zd.Zu 

V^Zu.Zd 



\Zu? + \Zd? \Zu? 



\Zd\ 



\/2zd.Zu 



\Zd? \Zu? + \Zd? 



(A7) 



TABLE I: Values of the fj. parameters (GeV^). 



/^ll ^12 




^^22 ^^23 




^33 


-7.0655 X lO-^ 1.6359 x 10^ 


-2.0184 X 10^ -1.2587 x 10"^ -1.7382 x 10^ 5.0687 x 10^ 


TABLE IL Values of the A parameters. 


All, 11 All, 12 Aiii3 


All,22 


All, 23 Aii_33 Ai2,12 


Al2,13 


Al2,21 


0.6385 -0.4227 -0.0347 


0.2500 


0.3128 0.8696 -0.0987 


0.2285 


0.3917 


Al2,22 Ai2,23 Ai2,31 


Al2,32 


Al2,33 Ai3,i3 Ai3,22 


Al3,23 


Al3,31 


-0.3132 0.1735 0.1780 


-0.0190 


0.4370 0.1852 -0.1830 


0.3268 


0.1620 


Al3,32 Ai3,33 A22,22 


^22,23 


A22,33 A23,23 A23,32 


A23,33 


A33,33 


-0.1566 -0.2230 0.2373 


-0.2803 


-0.1203 0.0536 0.4147 


0.1545 


0.5368 


CB one. The values of the vevs obtained are given below (we considered a CB minimum 
with the phase 6c equal to zero). With these vevs one finds the values of the potential quoted 


TABLE III: Values of the vevs 


for the normal and CB minima (GeV). 








Vcl Vc2 




a 


225.2135 -41.9564 


89.6355 


325.5199 343.9166 


17.9887 


in the main text. Using the methods developed in the first paper of Q| it is a simple matter 
to write down the squared mass matrices for the scalar fields. They are found to to have, 
other than the expected Goldstone bosons, only positive eigenvalues for both the CB and N 
stationary points, thus proving that both are minima. 
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